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Abstract. In this contribution we summarize recent results on tiie transport properties of strongly correlated dilute Fermi 
gases. We discuss the hydrodynamic equations in the normal phase and present new results on the structure of second order 
terms in the gradient expansion. We also discuss predictions from kinetic theory, and show how these predictions can be tested 
using experimental data on elliptic flow. We summarize current constraints on the shear viscosity to entropy density ratio r) /s. 
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In recent years a remarkable convergence has occurred between novel physical properties of the hottest matter, the 
quark gluon plasma, and the coldest matter, cold atomic gases, that can be studied in the laboratory [1]. Both types 
of matter, which are governed by almost scale invariant interactions, were found to behave like nearly perfect liquids 
whose shear viscosity is close to the limits thought to be set by quantum mechanics. Initial estimates showed that in 
both the quark gluon plasma and the dilute Fermi gas at unitarity the ratio of shear viscosity to entropy density satisfies 
rj/s <Q.5h/kB [2, 3, 4, 5]. This result is close to the value rj/s — h/{AnkB) which is found in the strong coupling 
limit of a large class of field theories that can be analyzed using the AdS/CFT correspondence [6, 7]. This suggests 
that strongly coupled quark gluon plasmas and dilute Fermi gases are more efficiently described in terms of a dual 
gravitational theory than in terms of a conventional quasi-particle model. 

In this contribution we summarize recent work on transport properties of dilute Fermi gases [8, 9]. These systems 
have been realized experimentally using optically trapped alkali atoms such as ^Li and ^"K. These two atoms are 
fermions because they posses a single valence electron and the nuclear spin is integer At very low temperature the 
atoms can be described as pointlike particles interacting via interatomic potentials which depend on the hyperfine 
quantum numbers of the valence electron. A Feshbach resonance arises if a molecular bound state in a "closed" 
hyperfine channel crosses the threshold of an energetically lower "open" channel. Because the magnetic moments of 
the states in the open and closed channel are in general different, Feshbach resonances can be tuned using an applied 
magnetic field. At resonance the two-body scattering length in the open channel diverges, and the cross section a is 
Umited only by unitarity, (y{k) ^4-7i/k^ where k is the relative momentum. 



At large distances and long times the dynamics of the dilute Fermi gas is described by hydrodynamics. For simplicity 
we will only consider the normal phase, corresponding to T > Tc ~ 0.147/?. Here, Tf = kp/{2m) is the Fermi 
temperature, kf — {37l^n)^^^ /{2m) is the Fermi momentum, and n is the density. In the normal phase there are five 
hydrodynamic variables, the mass density p = mn, the flow velocity v, and the entropy density s. These variables 
satisfy five hydrodynamic equations, the continuity equation, the Navier-Stokes equation, and the entropy production 
equation. 
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HYDRODYNAMICS 
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These equations close once we supply constitutive relations for the stress tensor n,y and the entropy current 7/ as well 
as an equation of state. The constitutive relations can be determined order by order in gradients of the hydrodynamic 
variables. The stress tensor is given by 

= pViVj + P5ij + SYlij , (4) 

where SYlij, the dissipative part, contains gradients of v,,p or s. At first order the dissipative contribution to the stress 
tensor is Slly = — TjCy with 

Oij = (ViVj + WjVi - ^5y(V,V,)) , (5) 

where t/ is the shear viscosity and we have used the fact that the bulk viscosity of the unitary Fermi gas is zero [10]. The 

entropy current is jf = VjS + 8 jf with 5 jf = Kqi and qi = V, log(r). Here, K is the thermal conductivity. Finally, the 
dissipative function R is given hy R = {r] /2)(y^ + Kq^. Scale invariance restricts the dependence of the shear viscosity 
and thermal conductivity on the density and the temperature. We can write 77 (n, 7") = an{y)n and K{n,T) = a„{y)n/m, 
where a„{y) and C7„(v) are functions of y = niT/n^^^. 

We have recently studied the constraints imposed by Galilean invariance and conformal symmetry on the structure 
of second order terms in the gradient expansion [11]. These terms are potentially important because the leading order 
result given in equ. (5), which imphes that the dissipative stresses are instantaneously determined by gradients of the 
velocity, can lead to acausal behavior [12]. We find that at second order 

+ Yi^iiTVj) T + yiV(iPVj)P + YiV^irVj^P + YA^(^j) T + 75^ {i^j)P , (6) 

where Xr, A,- and Ji are second order transport coefficients, Clij = V,v/ — Vyv; is the vorticity, and (.) denotes the 
symmetric traceless part of a second rank tensor. The most important part of this result is the structure of the terms 
proportional to Tr, which describe the relaxation of the dissipative stresses to the Navier-Stokes form. Conformal 
symmetry determines the particular linear combination of the comoving derivative of and (V*Vft)(7y that appear in 
this term. 



KINETIC THEORY 

Near Tc the transport coefficients t] and k are non-perturbative quantities that have to be extracted from experiment or 
computed in quantum Monte Carlo calculations. At high temperature (and at very low temperature, T <c Tc, see [13]) 
transport coefficients can be computed in kinetic theory. The shear viscosity was first computed in [14]. Here we will 
follow the recent work [15] which also considers the frequency dependence of the shear viscosity. In kinetic theory 
the stress tensor is given by 

where fp = f{t,x,p) is the distribution function of fermion quasi-particles and the factor 2 is the spin degeneracy. The 
equation of motion for fp is the Boltzmann equation. In order to extract the shear viscosity it is useful to consider the 
Boltzmann equation in a curved background metric gij. In this setting correlation functions of the stress tensor can be 
determined by computing variational derivatives with respect to the metric. The non-relativistic hmit of the Boltzmann 
equation in a curved space with metric gij is 

where F'^^. is the Christoffel symbol and C[f] is the collision integral. We consider small deviations from equilibrium 

and write / = fo + Sf with /o(p) = foip'p'gij/ {2mT)). We also write gij = 8ij + hij and linearize in hij and 5 /. We 
get 



This equation can be solved by going to Fourier space and making a suitable ansatz for df{(0,k;p). Once 5/ is 
determined we can compute Ily using equ. (7), and Tj(fi)) by matching the A: — )• limit to the hydrodynamic result 
equ. (4,5). The result is 



, , 7?(0) 15(mr)3/2 r 1 



(10) 



We observe that the shear viscosity is large in the weak coupUng Umit 0, and that t] decreases as the scattering 
length is increased. The shear viscosity saturates when the scattering length becomes comparable to the de BrogUe 
wave length a ^ Xds ^ (mr)^'/^. In this limit 7] only depends on Xds but not on the density or the scattering length. 
The frequency dependence of the shear viscosity is governed by the relaxation time x = {2>r])/(2e) [16]. We observe 
that relaxation is fast in the Umit where the viscosity is small. We also observe that the viscosity satisfies a sum rule 
which only depends on thermodynamic quantities, 

1 /"^ e 

- / dcoriio)) = -, (11) 

Jq 3 

where A ~ J is a cutoff that takes into account the breakdown of kinetic theory at very high frequency. The sum rule 
can be extended to A — >■ oo by subtracting the high frequency tail the shear viscosity, rj {(o) ^ rj {(o) — C/{15n^/rrm), 
where C is a temperature dependent constant known as the "contact" [17, 18, 19, 20]. 

There are several other interesting results that have recently been obtained in kinetic theory. Chao and Schafer 
determined the second order coefficients defined in equ. (6) [11]. They find, in particular, that Tr = (3Tj)/(2e) agrees 
with the relaxation time in equ. (10). Braby et al. determined the thermal conductivity [21], and Bruun as well as 
Sommer et al. computed the spin diffusion constant [22, 23]. The spin diffusion constant at unitarity is 

Ds = ^^i^] . (12) 



This result can be compared to the momentum diffusion constant — r]/p. Kinetic theory predicts that = 
457r^/^/(64v^m)(r/7»^/^, which implies that the ratio of Ds/Drj is independent of temperature. Near Tc, both 
diffusion constants are of order fi / m. A similar relation holds in the quark gluon plasma, where the ratio of the heavy 
quark diffusion constant Dq and the momentum diffusion constant D-q=r]/ {sT) is approximately constant [24]. 



ELLIPTIC FLOW 

The first experiment that demonstrated nearly perfect fluidity in the dilute Fermi gas was the observation of elliptic 
flow by O'Hara et al. [25]. The experiment involves releasing the Fermi gas from a deformed, cylindrically symmetric, 
trap. The density evolves as 

n{xux^,t) = , «o {xib^{t),xjbz{t)) , (13) 

where the transverse and longitudinal coordinate, /^x (0 = (0 ^^e scale factors, and «o {x± , ) is the equilib- 

rium density of the trapped system. The initial system is strongly deformed, Ar(P) = [{x^) / {x^)]^^^ «C 1. Hydrody- 
namic evolution converts the large transverse pressure gradient into transverse flow. As a consequence the aspect ratio 
ArIi) grows with time and eventually becomes larger than one, see Fig. 1. 

Viscosity slows down the transverse expansion of the system. In order to quantify the effect of shear viscosity we 
have to solve the Navier-Stokes equation for the expanding cloud [12, 26]. In general this has to be done numerically 
but for 77 = o^n with a„ ~ const very accurate semi-analytical scaling solutions can be found. For this purpose we 
make a linear ansatz for the force term fi = {ViP)/m = UiXi (no sum over i), where a,- = a,(?). The scale parameters a,- 
and bi are determined by the coupled equations 

bj_ 2/3 cOx fb±_ b. 



b^ , 4j3Aco^ / bi_ b 




FIGURE 1. The left panel shows data for the aspect ratio versus time, from [27]: Top Black, E = Q.bEf\ Pink, E = 2.3 Ef; Blue, 
E = 3.3Ef ; Bottom Red, E = 4.6Ef. Solid curves: Hydrodynamic theory with the viscosity as the fit parameter. The panel on the 
right shows Tj /s as a function of E/Ef. The red data points come from elliptic flow experiments, the blue data point were obtained 
from the damping of collective modes [27]. 



d, 

where OxjOj are the oscillator frequencies of the harmonic confinement potential (before the gas is released). The 
parameter j3 is defined by 

^ = (3A^aV/^ (£o/£f) ' ^^^^ 
where is the number of atoms, A = Ar{Q) the initial aspect ratio, and Eq/Ef the initial energy in units of 
Ef = {3NX)^^^Nco±. The initial conditions are b±{0) = b,{0) = 1, bj_{0) = k{Q) = 0, and aj^{0) = col, ^<i'^) = ^z- 
Dissipative effects fall into two categories. The terms proportional to j8 in equ. (14,15) correspond to friction - shear 
viscosity slows down the expansion in the transverse direction. The dissipative terms in equ. (16,17) describe reheating 
- shear viscosity converts some kinetic energy to heat which increases the pressure and eventually re-accelerates the 
system. Friction slows down the growth of (f ) as compared to the ideal evolution. Reheating reduces this effect by 
about a factor of 2. 

In order to make comparisons with data we have to take into account that a„ is not a constant. In [12, 26] we argued 
that a„ in equ. (18) should be interpreted as the trap average of the local ratio of shear viscosity over density. 




In the dilute corona of the cloud the local shear viscosity is independent of the density and equ. (19) is not well 
defined. This problem can be addressed by taking into account that the viscous stresses relax to the Navier-Stokes 
value on a time scale T that becomes large as the density goes to zero, see equ. (10). A relaxation model for (a„) was 
studied in [12]. An even simpler model can be constructed based on the assumption that the shear viscosity relaxes 
to its equilibrium value at the center of the trap and is proportional to the density in the dilute corona. This implies 
that ri{x) — Tj(0)(n(jc)/«(0)). This parametrization agrees with the relaxation model at the 30% level. It was used by 
Cao et al. [27] to analyze the data shown in Fig. 1. The hydrodynamic curves shown in Fig. 1 were obtained with 
rj — 77o(mr)^/^ and Tjo = 0.33. This agrees quite well with the prediction of kinetic theory rjo = 15/(32v^) ~ 0.26, 
see equ. (10). Once we have verified the T^^^ scaling behavior of the viscosity in the high temperature limit we can 
determine the shear viscosity to entropy density ratio all temperatures, see the right panel in Fig. 1. At very low 
temperature, T <Tf, dissipative effects in the elliptic flow experiment are very small and it is easier to determine the 
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viscosity from the damping of collective modes, see the blue points in Fig. 1. We observe that the minimum of tj/j 
occurs near the lowest temperature studied, and that (77 /s)min — 0.4. We emphasize that these values are averaged over 
the trap, and that more detailed studies are needed to find the true minimum of 77 /s. 

CONCLUSIONS AND OUTLOOK 

There are a number of issues that need to be addressed before an accurate value of r\ /s with fully controlled errors can 
be obtained. The most important of these is a better description of the transition from nearly perfect fluid dynamics in 
the center of the cloud to kinetic behavior in the dilute corona [28]. This problem can be addressed using the second 
order hydrodynamic equations discussed in this contribution. We would also like to develop tools that will allow 
us to perform calculations of transport properties in the strongly coupled regime. An important constraint for these 
calculations is provided by the sum rules presented in equ. (11). There have also been some very interesting attempts 
at extending the AdS/CFT correspondence to non-relativistic conformaUy invariant theories, see [29, 30]. Finally, we 
would like to understand whether nearly perfect fluidity in the dilute Fermi gas in the strongly interacting regime can 
be understood in terms of quasi-particle degrees of freedom. This question can be addressed by comparing different 
transport coefficients, Uke the ratio Dj/frj mentioned above, or by computing the spectral functions related to shear 
viscosity and diffusion. 
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